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Ricci curvature and the mechanics of soHds 

Michael Eastwoo 

Abstract 

We discuss some differential geometry pertaining to continuum 
meclianics and tlie route recently taken by D.N. Arnold, R.S. Falk, 
and R. Winther in deriving new improved finite element schemes in 
linear elasticity from constructions in projective geometry. 

Some vector analysis 

O 

CN ■ We start with some basics from vector analysis [3]. Let us write 

c^! di = d/dxi 82 = d/dx2 d^ = d/dx^ 

for the partial derivatives in M^. The gradient of a smooth function / defined 
(S| '■ on f/°P<=° C M3 ig the vector field 

grad/ = {dif,d2f,d3f) 

^ ' on ?7. If X = (Xi, X2, X3) is a smooth vector field on U, then 

curlX = (^2X3 - ^3X2, 93X1 - ^1X3, 9iX2 - ^2X1). 

It is readily verified that curl o grad = 0. Indeed, if f/ is a sufficiently simple 
set, such as a ball, then 

' fl) 



X = grad/, for some / <^=^ curlX = 0. 



^ . In effect, the curl of a vector field is the skew part of the 3x3 matrix (diXj) 

of partial derivatives. Let us instead consider the symmetric part 

^1 S = (S,,) = (|[9.X,+9,X,]) 

^J ', and ask for conditions that a given symmetric tensor field S = (Sjj) be of 

Q I this form. The answer is that S should satisfy the Saint- Venant equations 

curl curls = 0, where curl curl E is the symmetric matrix obtained by 



firstly regarding S as a row vector (whose entries just happen to be 



;_( ■ column vectors) to form curlS, 

• then regarding curl S as a column vector (whose entries just happen to 
be row vectors) to form curl(curlE). 

The statement ([1]) has a useful counterpart as follows. 



Sij- = lidiXj + djXi], for some X ^^ curl curl S = 0. (2) 



Indeed, we shall see that ([2]) can be deduced from ([T]). This, in turn, has 
consequences in the design of finite element schemes concerned with elasticity. 
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Ricci curvature in three dimensions 

Readers unfamiliar with differential geometry might omit this section on first 
reading. There is also a close link between ([2]) and Ricci curvature in three 
dimensions. Using the Einstein summation convention, if gij is a Riemannian 
metric with inverse g^^ , then the Ricci tensor Rij is the symmetric tensor 
given by 



dkT 



«j 



•^iJ- jk ' '- ij ^ mk 



p mp 

■ J- ifc J- jr, 



where r^'^ 



Ig^^idiQji + djgu-digij]. 



In three dimensions Rij = if and only if g^j is flat, meaning that there is 
a local change of coordinates that transforms g^ as a tensor into the flat 
metric 6ij (more specifically, there is a change of coordinates with Jacobian 
matrix J such that (gij) = J* J). The infinitesimal version of this statement 
is essentially (|2]). More precisely, if Ejj is an arbitrary symmetric tensor on 
^open ^ ]^3 ^T^^ ^g consider the metric gl, = Sij + eSjj where e is sufficiently 
small that gij is positive definite, then 



(curl curls) 



u 






where Gij is the Einstein tensor Rgij — 2Rij for R = g R^i- The Einstein 
tensor carries the same information as the Ricci tensor but has the advantage 
that the Bianchi identity simply says that Gij is divergence-free V^Gij = 0. 
Sure enough, one can readily verify that 5*(curlcurlE)jj = 0. 



Translation from (fTl) to (E 



Firstly, some convenient notation in three dimensions. Let us write eijk for the 
totally skew tensor with €123 = 1. It allows us to write (curlX)j = ei^^djX^. 
Now consider a pair F = (X^, Y^) of vector fields on U°^^^ C M^, regarded 
as a function with values in the vector space W = R^ © M^. If we define the 
gradient of F by 



grad 






djXi — eji^Yjr. 



d,Y, 



(3) 



and use this definition naively to compute the curl of a vector field with 
values in W, then we obtain 



curl 






Q Ojl^kt — Q ej£ ^km 



£«■' djTjke — S^j + SieEr 



ei^^djZu 



It is readily verified that curl o grad = 0. This says precisely that ([3]) defines a 
fiat connection, which enables one to deduce that if [/ is a sufficiently simple 
set, such as a ball, then 



^ = grad F, for some F 



curl\E' = 0. 



(4) 



To deduce ([2]), let us suppose that Sjj is symmetric and set 



^ 









SO that curl \1/ 







(curl curl E) 



if. 



If curl curls = 0, we immediately infer the existence of vector fields Xi and 
Y^ on t/ such that 



grad 



X, 



^ i.e. 



djXi — ejl^Yrr, 



d,Y, 



_ imf) V 



In particular, E^^ = |[9jX£ + S^Xj], as required. 

Continuum mechanics 

Although different words are used, Riemannian differential geometry in three 
dimensions is exactly what is needed to set up the mechanics of solids |1]. The 
metric tensor is known as the strain in continuum mechanics. The Einstein 
tensor is known as the stress. The Bianchi identity says that the stress tensor 
is divergence- free, interpreted as a conservation law in mechanics. Linearising 
around the flat metric gives the following complex of tensors on M^ 



displacement — )■ strain — t- 



Sij 



stress 



—7- load . 



fcm, £n 



(5) 



e/"aaE 



kiJi^mn 



where the displacement and load are vector fields whilst the stress and strain 
are symmetric 2-tensors. 

Finite element schemes 



We have already seen that the fiat connection ([3]) somehow embodies the 
operator S i— )■ curl curl S relating strain and stress in i^. More generally and 
precisely, the whole complex (^ may be derived from the connection (^. To 
do this, recall that the gradient operator (|3]) concerned functions with values 
in W = M^ © M^. Thus, we may write 



W 



grad 



curl 



®w 



div 



M3 
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R3 ©S^RS 
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R3©M3 
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R3 
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S2R3 © R3 


© 


M3 


R3©M3 


R3 



(6) 



where S^R^ denotes symmetric 3-tensors whilst skew 3-tensors are identified 
l^ using ej-'^. In this diagram, the spaces indicated thus I I are joined 



with 



by isomorphisms indicated thus /^. A simple diagram chase cancels these 



spaces and results in the linear elasticity complex (E]). In [I], Arnold, Falk, 
and Wintlier use a halfway-house complex 



obtained by cancelling only M M /^ R (g) M from (j6]), to construct 
new and stable finite element schemes for linear elasticity mimicking the 
previously known stable finite element schemes for the grad-curl-div complex. 



Projective geometry 

The connection ([3]) may be viewed as follows. Consider the unit three-sphere 
5''^ C M^. There is no difficulty in taking the gradient of a function F on 5*^ 
with values in the skew 2-tensors A^M^. However, each point on S^ is also a 
vector f G M^, which may be used to decompose these skew 2-tensors:- 

A^M^ = {u s.t. vAuj = 0}®{uj s.t. t; J w = 0} ^ M^ © M^ 

This decomposition does not see the sign of v and so descends to the quotient 
of S^ under antipodal identification, namely real projective 3-space MP3. 
The upshot is that the gradient of F may be written in terms of the intrinsic 
calculus on MP3 and viewed in a standard affine coordinate patch R^ "^ RP3. 
The result is ([3]). The construction of ([5]) from ([6]) is due to Calabi [2\. It 
may also be viewed as a geometric realisation of the Jantzen-Zuckerman 
translation principle from representation theory [6J and, as such, admits vast 
generalisation in the newly developed field of parabolic geometry [5] . 
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